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Abstract 

Neutrino-nucleon total cross sections for neutrino energies up to ultrahigh energies (UHE), 
= 10^^ GeV, are evaluated within the framework of the dynamical (radiative) parton 
model. The expected uncertainties of these predictions do not exceed the level of about 
20 % at the highest energies where contributions of parton distributions in the yet unmea- 
sured region around x ~ 10^* to 10^^ are non- negligible. This is far more accurate than 
estimated uncertainties of about 2^^ due to ad hoc extrapolations of parton distributions 
to a: < 10^^ required for calculating UHE cosmic neutrino event rates. 



Calculations of ultrahigh-energy (UHE) neutrino-nucleon total cross sections, relevant for 
neutrino astronomy, have been improved [jl], |^, ^, ^ by taking into account new high 
energy measurements of deep inelastic lepton-nucleon scattering (DIS) at DESY-HERA 
1^. These updated calculations were even further improved very recently not only 
due to additional small-x precision measurements at HERA [|^ which became available 
in the meantime, but mainly due to new ideas about the flux of neutrinos |§] from active 
galactic nuclei (AGNs) [§, gamma ray bursts (GRBs) or from decays of exotic heavy 
particles of generic top-down (TD) or topological defects models [|n|]. Most of these 
new developments refer to neutrino energies above about 10® GeV where (anti)neutrino- 
nucleon cross sections become sensitive to parton densities at ultrasmall values of Bjorken- 
X, X < 10~^, not accessible by present DIS experiments. Different assumptions about the 
X — s> behavior then lead to different cross sections. Thus a significant uncertainty 
in estimating the detectability of extraterrestial UHE neutrinos with present and future 
neutrino telescopes is due to the small-x extrapolations and the resulting uncertainty 
reaches typically a factor 2^^ around 10^° eV p|, |^. 

The highest presently available ep energy at HERA, ^Vs^ = 314 GeV, is rather small as 



compared to the envisaged UHE neutrino nucleon collision energies of up to about y/s^ = 
10® GeV. Estimates of the corresponding total '^v N cross sections thus afford either 
extensive, possibly unrehable, extrapolations [|l], ||, |^, ||, || of existing data (x > 10^^) and 
their respective fits, or the application of QCD inspired models which proved to provide 
reliable high energy predictions |T^, [13], in the past 0, 0. In the present work we shall 
adopt this second option and base our predictions on calculations within the framework 
of the radiative parton model |12|, |T3|, |l^, |15| which allows to calculate the small-x {x < 



10^^) behavior of parton densities from first principles, i.e. QCD dynamics, independently 
of any free (fit) parameters in the small-x region. These unique dynamical predictions 
result from valence-like gluon and sea input densities xg{x,Ql) ~ x"", xq{x,Ql) ~ x^' 
with a, a' > as X at some low momentum scale Qq ^ 0.5 — 0.6 GeV. The 
resulting small-x behavior is furthermore perturbatively stable and unique at the relevant 



momentum scales ~ M^. Here we shall apply the parton distributions in |T^, and in 



1 



particular the most recent ones in [IIH] which take also into account all recent high precision 



measurements at HERA g 0, to evaluate for 10^ < E(„) < lO^^ Q^y ^t 



highest neutrino energies, this requires the knowledge of parton densities down to a; ^ 10~^ 
(x > Myy/2MEi, due to the VT-propagator, with M being the nucleon mass) at scales 
~ M^. In ref. we have explicitly demonstrated that our next-to- leading order 



(NLO) results in this extremely small- a; region are perturbatively stable when compared 
with the leading order (LO) ones. It should be kept in mind that this stability refers 
always to measurable quantities like structure functions or cross sections rather than to the 



auxiliary, not directly measurable, NLO parton distributions ||T4|, [T^ since perturbative 
stability usually requires Q"^) 7^ f^'^{x, Q^) for f = q,q, Q- Furthermore, the NLO 

predictions are, in contrast to the LO ones, obviously rather insensitive to the specific 
choice for the renormalization scale appearing in as(/i|.) and for the factorization scale 
fip appearing in the parton densities f{x,fi'jp) [0. We shall therefore mainly concentrate 
on NLO analyses for predicting the various UHE total cross sections relevant for neutrino 
astronomy. 

The above results lend sufficient confidence in the reliability of the perturbatively 
calculated UHE N total cross sections, N = {p + n)/2, evaluated according to 

^''^'^(s) = /'dx f dy d^a''^'' /dxdy (1) 

Jo Jo 

with 

d^n'^^^ (72 (_) (_) (_) 

= ^i^ + ^ys/M^)-'[ii - y)F,^ + y'xF,^ ± y{i - l)xF,^ ] (2) 

where Fi = Fi{x,Q'^ = xys), s = 2ME(_) and Gp = 1.1663 x 10"^ GeV'^. It should be 
noted that the results for the total cross sections are insensitive to the lower integration 
limits in (1) since the W^-propagator in (2) restricts = xys to values around M^. 

fo, 

parton densities at their given input scale = Qq- The V N structure functions Fi in 
(2) can be decomposed as 

F, = Ff'""' + Fl' (3) 



Anyway , as soon as the integrations correspond to scales < Q^, we freeze the respective 



where F^^^^^^ refers to all light {u, d, s) contributions and Fj^ denotes all relevant contribu- 
tions related to heavy quark (c, b, t) transitions. 

For charged current (CC) uiN l^X reactions we have in LO-QCD 



Ff''''' = + d) + ^{d + u)\V^,\' + s\V^,\' 

pv,light _ 2rj,p'^^^^9ht 

pU,UgM ^ _^^^d) + {d + u)\V^d? + 2s\V^s? 

where u = u{x^Q'^) etc., except stated otherwise, s = s and the CKM mixing matrix 
elements are taken from |17]. The heavy flavor component F^ in LO consists of the 



charm contribution Ff according to the W^s — > c and W^d c transitions as described, 
for example, in |18[, and of the top-bottom contribution F/^ according to the W^g th 
fusion subprocess |jl9|, using a factorization and renormalization scale equal to {mt + mh) 



with rrit = 175 GeV and rrih = 4.5 GeV. [Note that the sign convention for F3 in [|T£ 
is opposite to the one in eq. (2)]. We shall compare these heavy quark contributions 
also with the less adequate, but frequently used (e.g., ^ |^, ^, massless approach 
(the so called 'variable flavor' scheme) where intrinsic 'heavy' quark densities are purely 
radiatively generated using the ordinary massless evolution equations, starting at Q = rrih- 
In this case, Fl^{x, Q^) in (3) is given by 

F^'^ = ^{d + u)\V,d\^ + s\V,s\'' + c + b{^,Q^ + m^,)\Vtb\'' 

F^''' = x{d + u)\Vcd\^ + 2xs\Vcs\^ + 2xc + 2^b{^,Q^ + m^t)\Vtb\^ (5) 



F^'"" = {d + u)\Vj' + 2s\Vj'-2c + 2b{^,Q' + mi)\Vt 



tb\ 



with ^ = x{l + mf/Q'^) and where we neglected ml/Q'^ contributions as well as t{x, Q"^). 
We choose the scale + mf in b, since it resulted in the best perturbative LO/NLO 
stability in the cs sector [ffH]. For this latter analysis in the so called 'variable flavor' 



scheme we shall use the dynamical parton densities of ref. [0. The relevant structure 
functions for CC uN reactions in (2) are simply obtained from the above expressions via 
Fr2 = FUq ^ q) and F^ = -FI^{q ^ q). 



3 



In NLO-QCD all above LO parton densities have to be replaced by the NLO ones [|I| 



15 1 , with the appropriate addition of convolutions with the fermionic Wilson coefficient 



Cq and the NLO contribution Cg®g, as summarized for example in [20|. It turns out that 



these additional convolutions contribute less than about 2% and thus these complications 
can be safely neglected for NLO calculations of high-energy total N cross sections. 
The heavy flavor piece receives additional massive 0{as) contributions according to 
the W^g — > cs, W^s' gc, etc., subprocesses |2l], ^ which we include according to 
the (MS) results of ||TE|. The fact that the tb production has so far only been calculated 
in LO {W^g tb) is of minor importance, since the cs sector dominates over the much 
heavier tb one. It should be noted that this more appropriate way of calculating the heavy 
quark contributions in fixed order perturbation theory without neglecting mass effects. 



usually referred to as 'fixed flavor' factorization scheme (also adopted in 0)' results 



in perturbatively stable LO/NLO QCD-predictions |jT8[ even for the highest energies 



attainable with present cosmic ray shower experiments {^/s = 10^ — 10^ GeV). 

For completeness we also summarize briefly the structure functions in (2) for neutral 
current (NC) reactions, where obviously Myy — Mz- For uN vX the light quark 
contributions in LO are 

2Fr'^^'* = \{u + u + d + d)iy: + Al) + \{u + u + d + d + As)iy! + Al) 

pu^ght ^ 2xF'('^''^' (6) 
^p.,ught ^ 2{u, + d,){V^A^ + VdA^) 

where g„ = qy{x,Q'^) = q - q and K = | - |sin^6'vy, Vd = -| + |sin^6'^y, A^ = 
—Ad = \ with sin^ dw = 0.232. The (massive) heavy flavor contributions to F^'2 derive 
from the subprocesses Z^g hh, h = c,b and are given in m, with F.^'" = 0. (The 
small contribution from ti production will be neglected). If heavy quarks are treated less 
adequately as massless intrinsic partons, we have instead of eq. (5), 

2F,^'' = 2c (K' + Al) + 2b {V! + Al) 

F!^'^ = 2xFl'^ , F3"''' = . (7) 



In this massless 

Q2) etc., we use again the 'variable flavor' parton densities 



of ref. |13|. For NC uN vX reactions, one simply has = F^. In NLO all above LO 



parton distributions have to be replaced by the NLO ones |14, 15 1, with the additional 



convolutions Cq® q and Cg® g |]20[ being again negligible as for CC reactions.[| 



Our relevant most recent NLO (MS) sea and gluon densities [|I^ at = are 
shown in fig. 1, where they are also compared with our previous GRV 94 [1^ predictions^] 
which are slightly steeper for x < 10"^, relevant for uN cross sections at E,, > 10^ GeV. 
For illustration we also show the input of the valence-like NLO GRV 98 densities at the 
input scale Qq = 0.40 GeV^ which becomes vanishingly small at x < 10~^, particularly 
for the gluon input. This illustrates the purely dynamical (i.e. parameter-free) origin of 
the very small- a; structure of sea quark densities and in particular of the gluon density 
at > Qq which dominates the evolution of q{x, Q"^) in the small-x region. Also 
noteworthy is the stability oiu + d dX Q"^ = 3> Ql with respest to our previous GRV 
94 results. For comparison the expectations based on the CTEQ3-DIS and CTEQ4- 



DIS parton densities]^ are shown as well. Whereas the CTEQ3 extrapolations seem 
to overestimate the very small-x region, x < 10~^, the CTEQ4 ones are in reasonable 
agreement with our dynamical expectations. 

The resulting CC total uN cross sections are presented in fig. 2 which demonstrates 
the stability of the QCD predicted cross sections in the UHE region. The heavy quark 
contributions are also shown separately: The cs' sector contributes significantly, whereas 



^ Since the fully massive NLO contributions to hh have not yet been calculated, we use the 

LO(as) expressions for the NLO analysis as well. (The factorization scale is taken to be = Ami 
which provides the best perturbative LO/NLO stability for photon- induced processes 7*3 hh, etc. [g3| 
- a stability which holds down to x — 10^^). The same holds for the above CC contribution due to 
W^g —^ tb. This approximation is of minor importance, since the main contributions come, in both CC 
and NC cases, from the perturbatively stable light u, d, s quarks and the cs sector. 

^ The simple analytic parametrizations of the LO and NLO QCD-cvolution GRV 94 predictions for 
parton densities presented in [ p^ are sufhciently accurate down to the lowest relevant x ~ 10~*. Moreover, 
the NLO(DIS) results are practically the same as the NLO (MS) ones for GRV 94 as well as for the 
more recent GRV 98 distributions. 

These densities have been extrapolated to low-x (x < 10~^) beyond the validity of the fitted 
parametrizations {x > 10~^) using some reasonable analytic functions. As in the case of the CTEQ3 



densities |24 , the CTEQ collaboration provided us with similar extrapolation functions for the CTEQ4 



distributions p5|. We thank Wu-Ki Tung for a helpful correspondence. 
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the much heavier tb contribution is small. These exactly calculated massive results are also 
compared with the results obtained from the less adequate description of treating heavy 
quarks as massless intrinsic partons according to eq. (5), using the GRV 92 densities |[T3| . 
Although this latter simplified massless approach to 'heavy' quark effects ot'enestimates 
their contributions, their difference with respect to the exact massive treatment becomes 
rather marginal at scales as large as ^ M^. It should be mentioned that in practically 
all conventional approaches where parton densities are obtained from fitting appropriate 
parametrizations to presently available DIS data, to be discussed below, heavy quarks are 
treated as massless intrinsic partons. This is also the case for the GQRS 98 results 
shown in fig. 2 which are based on the CTEQ4-DIS parton distributions as the nominal 
distributions for calculating neutrino-nucleon cross sections and estimating astrophysical 
UHE neutrino event rates. It should be kept in mind that these results have been obtained 
by using some ad hoc 'brute force' extrapolation to x < 10~^, i.e. by extrapolating the 
CTEQ4-DIS densities below x = 10^^ with the same power in x predicted at a; = 10~^ 
and = . This overestimates the light sea quark densities at x ~ 10~® by 

almost 20% as compared to the CTEQ extrapolations shown in fig. 1 and thus also 
overestimates the cross section by the same amount at Ei, ~ 10^^ GeV. Nevertheless, the 
expected total cross sections agree, somewhat accidentally, rather well with our dynamical 
predictions and are thus perfectly legitimate for practical applications. 

To test the amount of extrapolation in x and involved in these calculations of total 

cross sections, fig. 3 shows the contribution to cr^^ from different regions of x, i.e. by 

choosing a finite lower limit of x-integration x > Xmin in eq- (1)^ Even at highest neutrino 

energies the VT-propagator in (2) constrains the relevant values of x to be larger 

than about 10~^. For similar reasons the relevant is restricted to values near for 

E,y > 1 TeV. Since our dynamical QCD small- a; predictions for x < 10^^ agree with all 

present DESY-HERA measurements down to a; ~ 10~^, as discussed at the beginning, 

our dynamical predictions for UHE z/(i/) nucleon cross sections, which are dominated by 

X > 10~^ — 10~^, appear to be reasonably reliable. 

* We thank Ina Sarcevic for a clarifying correspondence on this point. 



6 



Similar results hold for CC vN and for the NC v{i')N total cross sections. Our 
final dynamical NLO high energy (small-x) predictions, based on our GRV 98 parton 
distributions [|l^, are given within 5 to 10% by 
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In fig. 4 we compare our most recent predictions for representative example, 

based on the GRV 98 distributions, with calculations based on different sets of parton 
densities and small-x extrapolations. It should be again pointed out that all conventionally 
fitted sets of parton densities are extrapolated to low-x (x < 10^^) beyond the validity of 
the fitted parametrizations at x > 10"^ using either some 'reasonable' analytic functions^ 
2j, 1^ or, in most ixed slope of xq{x, M^) at the lowest x permitted by present 



HERA experiments {x ^ 10~^). As already demonstrated in fig. 1, the dynamical GRV 
94 densities result in very similar high energy predictions as the GRV 98 ones, with GRV 
94 cross sections being about 5 to 10% larger for E^, > 10^ GeV. This is due to the fact 
that the GRV 94 sea and gluon densities are slightly steeper at small-x than the GRV 
98 ones. Notice that such a difference is within the typical uncertainty of about 20% at 
X ~ 10~8 due to factorization scale ambiguities The same holds true for the results 
in fig. 4 based on the CTEQ4-DIS densities as used and extrapolated by GQRS 98 



and the CTEQ3-DIS distributions |2J] which provide equally reasonable and acceptable 
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estimates of total neutrino cross sections 0, except perhaps at highest energies where 
GQRS 98(CTEQ4-DIS) possibly slightly underestimates them. The situation changes if 
we turn to the results based on the remaining extrapolated parton parametrizations in fig. 
4: The FMR extrapolations and the fixed-power extrapolated MRS-G distributions 
p7| result in noticeably larger cross sections, up to almost a factor of 2, at large energies as 
compared to our dynamical QCD extrapolations and predictions. On the other hand, the 
EHLQ-DLA distributions [^, extrapolated to x < 10~^ using the 'double logarithmic 
approximation' (DLA), strongly underestimate the cross sections almost throughout the 
whole high-energy range shown in fig. 4. 

These results indicate that UHE (anti)neutrino nucleon total cross sections can be 
calculated with an uncertainty of about ±20% at highest neutrino energies of 10^^ GeV 
which requires a reliable knowledge of parton distributions at a; = 10~® to 10~^ and 

= My^ [0. This is considerably more accurate than the estimated uncertainty factor 
of 2^^ derived from ad hoc extrapolations of fitted parton distributions to very small- 
X beyond the region of their validity [|], ^. Thus all recent estimates of cosmic UHE 
neutrino event rates |Q, |^, ^, |^ based on the dynamical GRV densities |T3| , |T^, |T^] and the 



properly extrapolated CTEQ distributions or parton densities with an effective 

similar small-x behavior, appear to be realistic and accurate to within 20%; this is in 
contrast to estimates based on extrapolated versions of FMR [|l| and MRS-G [|2|, ^ or 
MRS-D'_ P] distributions, for example, which result in too large rates at highest neutrino 
energies of about 10^^ eV. These conclusions hold of course mainly for downward event 
rates, since the upward muon event rates are obviously rather insensitive to the particular 
choice of parton distributions due to the compensating attenuation of neutrinos as they 
pass through the Earth |^ . 
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Figure Captions 



Fig. 1 The light quark sea and gluon distributions in NLO a.t Q = M^r. The dynamical 



GRV 94 [14] and GRV 98 [15| predictions have been obtained from evolving the 
respective valence-like input at = Ql which is shown for illustration for the 
GRV 98 densities by the two curves at Ql = 0.40 GeV^. For comparison, the 



extrapolations of the CTEQ3 and CTEQ4 [^'^ distributions are shown as well. 

Fig. 2 Charged current uN cross sections as calculated in NLO. In the case of GRV 98 
densities heavy quark contributions are treated in a fully massive way based 



on fixed-order perturbation theory as explained in the text, whereas the GRV 92 
distributions |13| refer to a massless (resummed) treatment of heavy quark flavors 
according to eq. (5). The GQRS 98 calculation is based on the CTEQ4-DIS 
densities extrapolated to x < 10~^ using a flxed power in x as given at x = 10~^ 
and = M^; this extrapolation method results in cross sections about 20 % larger 
at = 10^2 GeV than the CTEQ4-DIS densities shown in flg. 1. The HERA 
measurement is taken from ref. [^. 

Fig. 3 The dependence of cx^^ on the small-x region, i.e. on the lower integration limit 



Xmin < X in eq. (1), using the NLO GRV 98 parton distributions [|T^] 



Fig. 4 Ratio of charged current cross sections using different parton distribution functions 
(PDF) as compared to the GRV 98 distributions [l|]. The PDF refer to GRV 94 
H, GQRS 98 (CTEQ4-DIS) GQRS 96 (CTEQ3-DIS, MRS-G) @, FMR 
and EHLQ-DLA §. 
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